Thermodynamics of the spin-flop transition in a quantum XY Z chain 
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A special limit of an antiferromagnetic XYZ chain was recently shown to exhibit interesting bulk 
as well as surface spin-flop transitions at T = 0. Here we provide a complete calculation of the ther- 
modynamics of the bulk transition using a transfer-matrix-renormalization-group (TMRG) method 
that addresses directly the thermodynamic limit of quantum spin chains. We also shed some light 
on certain spinwave anomalies at low temperature predicted earlier by Johnson and Bonner. 
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There has been a revival of interest in bulk and sur- 
face spin-flop transitions following, some recent experi- 
mental work on Fe/Cr multilayers .El These are effectively 
described by classical spin chains characterized by anti- 
ferromagnetic exchange interaction in addition to single 
ion anisotropy. It is then natural to raise similar ques- 
tions in the context of quantum spin chains which are 
more appropriate for the study of quasi-one-dimensional 
crystalline magnetic systems. p 

Indeed, in a recent communication,^ this issue was 
studied within a special limit of the spin— i XYZ chain 
defined by the Hamiltonian 
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where A > 1 and the operators Tf satisfy the standard 
spin commutation relations; hence (|^) may be thought of 
as the Hamiltonian of a ferromagnetic XX Z chain in a 
staggered magnetic field. A more physical interpretation 
is obtained by the canonical transformation = T^, 
Sy = (-l)^r/, 51 = (-l)^r/ which reduces (0) to a 
special limit of an antiferromagnetic XYZ Hamiltonian 
in a uniform field H . In view of the above dual interpre- 
tation we consider in parallel the two special operators 



(2) 



e=i 



The operator r is not endowed with a simple physical 
meaning but commutes with Hamiltonian dl^) and thus 
provides a very useful classification of states. In contrast, 
the operator AI does not commute with the Hamiltonian 
but represents the physical magnetization, a quantity of 
special interest in the following. 

Although the main objective of this paper is to study 
the thermodynamic limit A — > oo, some issues become 
clear by approaching that limit through a finite periodic 
chain with an even number of sites A — 2N. The eigen- 
values of T are then given by t = 0, ±1, . . . , ±7V and split 



the Hilbert space into 2N + 1 sectors. The two extremal 
sectors r = ±N contain only one state each, which is an 
exact eigenstate of the Hamiltonian with energy E = Q 
for any strength of the applied field H. In the original 
spin language these are the two completely polarized Neel 
states. To study their stability at finite H we also con- 
sider one-magnon excitations, with r = — 7V-|- 1 or TV — 1, 
whose energy eigenvalues are the same for both sectors 
and are given explicitly by 



Ek^A± Vcos2(fc/2) + 7f2, 



G± = A ± Vl + H^, 
(3) 



where A: is a sublattice crystal momentum. In Eq. (|3|) we 
also display the energies of the k — modes, denoted by 
G±, which will be referred to as the magnon gaps and 
are depicted in Fig. 1 as functions of the applied field. 
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FIG. 1. Field dependence of the magnon gaps for the 
speciflc anisotropy A = 3/2 for which the critical field is 
Hh = 1.118. The skew dashed line represents the Ising asymp- 
tote H - A. 

It is clear that the lowest gap closes (G- = 0) at the 
critical field 

Hb - VA2-1 (4) 

beyond which the Neel states are no longer the lowest- 
energy states and the system undergoes a bulk spin-fiop 
(BSF) transition. 
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The nature of this T = phase transition is actually 
more interesting than indicated by the preceding argu- 
ment. In fact, the lowest-energy states of all sectors be- 
come degenerate at the critical field, with energy E — 0, 
and the coixesponding eigenstates can be constructed 
analyticallyE For H > Ht, the r = sector prevails in 
the sense that it contains the unique absolute ground 
state. Accordingly the first excited states are the lowest- 
energy states of the r = ±1 sectors and are degenerate. 
The corresponding magnon gap, denoted by G in Fig. 1, 
was computed via a Lanczos algorithm, on finite peri- 
odic chains with A < 22, complemented by straightfor- 
ward Richardson extrapolationE The detailed numerical 
results indicate that the gap G might vanish through an 
essential singularity at Hb in the thermodynamic limit. 
Extrapolation becomes completely unnecessary for fields 
in the region H }^ A where the gap G approaches the 
Ising asymptote iJ — A. One would expect that the low- 
est gaps and G dominate the low-temperature ther- 
modynamics, in the respective field ranges, an issue that 
turned out to be more intricate than normally antici- 
pated. 

In order to prepare the discussion of thermodynamics it 
is also useful to calculate the magnetization Af at T = 0. 
The magnetization vanishes for H < Hb but exhibits a 
finite jump at the critical field which can be calculated 
analytically. For H = Hb the expected value of M in the 
ground state \iPt) of each sector r is 
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which generalizes the r = result quoted in Ref. 2. Fog. 
any fixed r a simple application of the Laplace methodu 
yields the asymptotic expansion 
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which can be used to extract the thermodynamic limit. 
Since the r = sector contains the absolute ground state 
just above Hb, the magnetization jump at the critical field 
is given by 
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A more subtle quantity is the T = average magnetiza- 
tion at the critical point calculated from 



fib = hm . „ M. 

^ ' T=-N 



(8) 



The asymptotic expansion cannot be employed in Eq. 
(0) because the latter contains terms with values of r that 



are comparable to N. Nevertheless the explicit result (|^) 
may be inserted in Eq. (H) to numerically estimate fib at 
large N. 

For H > Hb the T — magnetization is not known 
analytically and we have again resorted to the Lanczos 
algorithm. At our maximum size A = 22 the Lanczos 
result for the magnetization jump at the critical field dif- 
fers from the analytical prediction (|^) by about 5%, a 
difference that is rectified by Richardson extrapolation 
to an accuracy about one part in a thousand. Hence we 
have applied the same extrapolation for H > Hb and 
the result is depicted by a solid line in Fig. 2. Again, 
extrapolation becomes unnecessary for H ^ A. 
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FIG. 2. Field dependence of the average magnetization per 
site fi = M/A for r = (solid line) and T = 0.1 (dotted line). 
The solid circles represent the values of the T — Q magneti- 
zation just below Hb (fi — 0), right at Hf, {fit, = 0.143), and 
just above Hb [fio ~ 0.224). The meaning of the open circles 
is discussed in the text. 

We thus arrive at the main point of this paper, 
namely the calculation of thermodynamics via a TMRG 
algorithmo which has alrjeady been applied to the study 
of quantum spin laddersO. One of the distinct features 
of the method is that it directly addresses the thermody- 
namic limit A — > oo. We shall not present here numerical 
details but merely discuss some important results. 

For instance, the temperature dependence of the mag- 
netization is shown in Fig. 3 for a number of field values. 
The magnetization vanishes for all temperatures at van- 
ishing field. For finite fields in the subcritical region, 
H < Hb, the magnetization again vanishes at T = 0, as 
expected, but develops a maximum at some finite tem- 
perature. Right at the critical field, H = Hb, the T = 
limit of the calculated curve is consistent with the value 
fib — 0.143 of Eq. (||), applied for A — 3/2, whereas just 
above Hb the T — limit is consistent with the value 
^0 = 0.224 of Eq. (0). For supercritical fields, H > Hb, 
the low-temperature limiting values of fi extracted from 
Fig. 3 are depicted by open circles in Fig. 2 and are 
thus seen to be in excellent agreement with our indepen- 
dent Lanczos calculation of the magnetization at T = 0. 
Similarly the results extracted from Fig. 3 at the spe- 
cific temperature T = 0.1 were used to calculate the field 
dependence of the magnetization at this temperature, a 
result that is shown by a dotted line in Fig. 2 and il- 
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lustrates the manner in which the T = magnetization 
jump at the critical point is smoothed out at finite tem- 
perature. 
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FIG. 3. Temperature dependence of the magnetization for 
various field values throughout the BSF transition. The solid 
circles on the jj, axis correspond to the same values as those 
of Fig. 2. 

We next turn our attention to the specific heat. In Fig. 
4 we compare the TMRG result at vanishing field with 
a finite-size calculation for chains with A = 10, 12, and 
14 for which a complete numerical diagonalization of the 
Hamiltonian is possible. This comparison is surprising in 
that the trend of the finite-size results does not seem to 
be consistent with the calculated thermodynamic limit. 
We have thus naturally questioned the validity of our 
TMRG calculation. However this special case was also 
considered in Fig. 4b of a paper by KliimperO whose 
numerical method is again based on a transfer matrix 
but relies heavily on the complete integrability of model 
(|l|) at vanishing staggered field. Direct correspondence 
with the above author established that our result agrees 
with his throughout the temperature range considered. 
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FIG. 4. Temperature dependence of the specific heat per 
site C at vanishing field [H — 0). 

It should be added here that the TMRG method does 
not rely on complete integrability and is thus more flex- 
ible; e.g., it can be applied for the calculation of the 
thermodynamics at any finite staggered field for which 
model (|l]) is not known to be completely integrable. 

The "anomalous scaling" observed in Fig. 4 for H ^ 
persists for nonvanishing fields throughout the BSF tran- 



sition but gradually disappears in the "no-scaling region" 
H ^ A where the correct thermodynamic limit is prac- 
tically reached by very short chains, as short as A = 4. 
In any case, the TMRG calculation of the temperature 
dependence of the specific heat is illustrated in Fig. 5 
for various field values. The main feature of this figure 
is that the specific heat develops a double peak for fields 
in the vicinity of the critical point Hi,. Furthermore the 
low-temperature behavior appears to be generally consis- 
tent with the field dependence of the magnon gaps shown 
in Fig. 1. 
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FIG. 5. Specific heat for various field values throughout the 
BSF transition. 

One would expect that the low-temperature specific 
heat is correctly predicted by a dilute-magnon or spin- 
wave approximation, a long cherished assumption in con- 
densed matter physics. To check this assumption we first 
consider the case of vanishing field for which our model 
is formally identical to the ferromagnetic XJf,Z chain 
extensively studied through the Bethe AnsatzO'Ll. At suf- 
ficiently low temperature the spinwave approximation of 
the specific heat should read 



C 



Gl exp(-G_/T) 
(27rr3)i/2 ' 



G_ = A - 1, 



(9) 



where G- is the lowest magnon gap at vanishing 
field. Equation (||) suggests considering the quantity 
— Tln(T'^/^C) which should interpolate linearly to the 
magnon gap G_ at T = 0. Yet a comparison of the spin- 
wave prediction with the TMRG calculation shown in 
the A = 3/2 entry of Fig. 6 reveals a sharp disagreement 
even at the lowest temperature accessible by our method. 
On the other hand, one can show that the spinwave ap- 
proximation agrees well with the finite-size results for 
A = 10, 12, 14 given in Fig. 4 restricted to the tempera- 
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ture range of Fig. 6. Clearly then the anomalous scaling 
noted earlier is intimately related to the apparent failure 
of the dilute-magnon approximation. 
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FIG. 6. Comparison of the calculated specific heat (solid 
line) at vanishing field {H — 0) with the spinwave approxi- 
mation (9) (dotted line) and the Johnson-Bonner prediction 
(10) (dashed line). 

In order to understand this situation we now invoke an 
asymptotic result obtained for the ferromagnetic XX Z 
chain by Johnson and BonnerQ who predict that the low- 
temperature specific heat is more appropriately described 

by 

G^exp(-G_/T) G?cxp(-Gi/r) 
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gion T <^ S is difficult to approach by the inherently 
finite-temperature TMRG algorithm. This explains the 
apparent failure of spinwave theory demonstrated in the 
A = 3/2 entry of Fig. 6. However, when both terms of 
Eq. (nw are included, the agreement with our TMRG 
result is obviously very good. The picture becomes more 
transparent in the A = 5/4 entry of Fig. 6 where the 
differential gap 6 = 0.125 is greater and thus the region 
T ^ 6 becomes accessible to TMRG, albeit somewhat 
marginally. Also interesting is the result for the critical 
anisotropy A = 5/3 shown in Fig. 6, where the failure 
of spinwave theory becomes complete, whereas our result 
continues to agree with the Johnson-Bonner prediction 
( p^ ) . Finally we have examined the case of a supercriti- 
cal anisotropy, A = 2, with a similar conclusion. 

At finite (staggered) field our model is not equivalent to 
the ferromagnetic XX Z chain and thus the finite-field re- 
sults of Ref. 7 are no longer applicable. We do not know 
at this point how to generalize Eq. (p^ ) to account for a 
staggered field, especially because complete integrability 
seems to be lost. Numerical investigation of this issue 
suggests that spinwave anomalies persist in the subcrit- 
ical region H < Hb while normal spinwave behavior is 
restored for H > Hh. In the latter region the quantity 
—Tln{T^^^C) interpolates linearly to the magnon gap G 
shown in Fig. 1. 

To summarize, we have presented a reasonably com- 
plete theoretical description of the thermodynamics of 
the spin-flop transition for Hamiltonian (|^). Our explicit 
results would be directly relevant for the analysis of ac- 
tual experiments, provided that a quasi-one-dimensional 
magnetic system is found that is described by our model 
Hamiltonian at least approximately.a Perhaps equally im- 
portant is the overall conclusion that the TMRG method 
proves to be reliable even under stringent conditions. 
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(10) 



where a new gap Gi is potentially important. This gap 
originates in bound multimagnon or domain- wall states, 
including the notorious factor 1/2 familiar from ear- 
lier disc-ussions of the antiferro mag netic Ising and XX Z 
chains.u The two terms in Eq. ( [10[ ) may then be referred 
to as the magnon and Ising contributions, respectively. 
The two gaps G_ and Gi become equal at the critical 
anisotropy A = 5/3 and are ordered as G_ < Gi or 
G_ > Gi for A < 5/3 or A > 5/3. 

Therefore, when 1 < A < 5/3, the magnon contri- 
bution in Eq. (|lO|) dominates for sufficiently low tem- 
perature, practically in the region T <C Gi — G_ = S. 
For A = 3/2 one finds that 5 — 0.06 and hence the re- 
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